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ABSTRACT 

We have used the Keck I Telescope to measure i\~-band surface brightness fluctuations (SBFs) of 
NGC 4874, the dominant elliptical galaxy in the Coma cluster. We use deep HST WFPC2 optical 
imaging to account for the contamination due to faint globular clusters and improved analysis techniques 
to derive measurements of the SBF apparent magnitude. Using a new SBF calibration which accounts 
for the dependence of iGband SBFs on the integrated color of the stellar population, we measure a 
distance modulus of 34.99 ± 0.21 mag (100 ± 10 Mpc) for the Coma cluster. The resulting value of 
the Hubble constant is 71 ± 8 km s" 1 Mpc" 1 , not including any systematic error in the HST Cepheid 
distance scale. 

Subject headings: distance scale — galaxies: distances and redshifts — galaxies: elliptical and 
lenticular, cD — galaxies: individual (NGC 4874) — infrared: galaxies 
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IR wavelengths. K-b&nd (2.2 /im) SBFs have been mea- 
sured for galaxies in several clusters (Luppino & Tonry 
1993; Pahre & Mould 1994; Jensen et al. 1998, 1999; Mei 
et al. 2001; Liu et al. 2001). 

The core of the Coma cluster is dominated in luminos- 
ity by two supergiants, the elliptical NGC 4889 (Bt = 
12.53 mag) and the cD NGC 4874 (B T = 12.63 mag; 
de Vaucouleurs et al. 1991). Like many Abell clusters, 
Coma possesses significant substructure (e.g., Davis & 
Mushotzky 1993; Mohr et al. 1993). However, NGC 4874 
appears to be special for two reasons: it is located at peak 
of the diffuse X-ray emission (White et al. 1993), and it 
has a strong nuclear radio source characteristic of many 
central giant ellipticals (Harris 1987). Therefore, it is con- 
sidered to reside at the heart of the gravitational potential 
of the cluster (Baier et al. 1990; Colless & Dunn 1996). 

In this paper, we measure the distance to NGC 4874 us- 
ing X-band SBFs. We have recently developed new theo- 
retical models for SBF studies (Liu ct al. 2000, hereinafter 
Paper I). Also, we have completed a new calibration of K- 
band SBFs using a large sample of early-type galaxies in 
nearby clusters (Liu et al. 2001, hereinafter Paper II). We 
have found that if-band SBFs vary considerably between 
galaxies but, like /-band SBFs, do so in a predictable fash- 
ion as they are correlated with V — I c galaxy color. In 
this paper, we use these results with new SBF data from 
the Keck I Telescope to derive an accurate distance to 
NGC 4874 and consequently a measurement of Hq. 

2. observations 

We observed NGC 4874 at the 10-meter Keck I telescope 
using the facility instrument NIRC (Matthews & Soifer 
1994) with the standard if-band filter (2.0-2.4 ^m). The 
camera uses a Santa Barbara Research Corporation 256 x 
256 InSb array and has a 38" x 38" field. Table 1 summa- 
rizes the observations. Both runs were photometric, with 

1 Based in part on observations obtained at the W.M. Keck Observatory, which is operated as a scientific partnership among the California 
Institute of Technology, the University of California, and the National Aeronautics and Space Administration. 

2 Currently Beatrice Watson Parrent Fellow, Institute for Astronomy, University of Hawai'i, 2680 Woodlawn Drive, Honolulu, HI 96822. 
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1. INTRODUCTION 

The Coma cluster (Abell 1656) is an important rung in 
the cosmic distance ladder. It is the nearest of the very 
rich Abell clusters (richness class of 2 [Abell 1958]) and 
thus contains many luminous galaxies which can be tar- 
geted for distance measurements. Also, the cluster redshift 
(«7000 km s" 1 ) is large enough that its cosmological re- 
cession is not significantly perturbed by smaller-scale pe- 
culiar velocities (e.g. Strauss & Willick 1995). For both 
these reasons, several methods have been used to measure 
the distance to Coma and thereby determine the Hub- 
ble constant {Hq). These include Type la supernovae 
(Capaccioli et al. 1990), the Tully-Fisher relation (Gio- 
vanelli et al. 1998), the Fundamental Plane (Jorgensen 
et al. 1995; D'Onofrio et al. 1997), the globular cluster 
luminosity function (e.g., Kavelaars et al. 2000), and sur- 
face brightness fluctuations (Thomsen et al. 1997; Jensen 
et al. 2001). 

Surface brightness fluctuations (SBFs) are an appeal- 
ing distance indicator because they have a well-understood 
physical basis — they arise from Poisson statistical fluctu- 
ations in the number of stars in a resolution element (i.e., 
the seeing disk). Over the past decade, optical SBFs have 
been used to measure distances to elliptical galaxies and 
the bulges of spiral galaxies (Jacoby et al. 1992; Blakeslee 
et al. 1999). /-band (0.8 /*m) SBFs vary between galax- 
ies by up to one magnitude, but the variations are well- 
correlated with V — I c color so they can be compensated. 
More recently, near-infrared (IR) SBFs have been shown 
to be useful distance indicators. Because cool giant stars 
dominate the spectral energy distributions of ellipticals, 
SBFs are brighter in the near-IR and hence can be de- 
tected at greater distances. Unresolved globular clusters 
are the dominant contaminant for SBF measurements at 
large distances, but since the clusters are bluer than the old 
stars in ellipticals, the contamination is much smaller at 
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clusters and background galaxies fainter than the detection 
limit are unresolved point sources in our images. There- 
fore, they contribute to the fluctuations of the galaxy sur- 
face brightness. The effect is substantial, because nearly 
the entire globular cluster population of the galaxy is un- 
resolved in the if-band images (see below). In order to 
avoid this contamination, we identified the positions of 
the globular clusters in our K-band data using archival 
deep HST U-band imaging of NGC 4874 (Kavelaars et al. 
2000), which reach down to the peak (V s=s 27.9 mag) of 
the globular cluster luminosity function (GCLF). 

We then constructed a software mask defining an an- 
nular region and excluding the globular clusters. For the 
1998 Keck data, we used two annuli, 3-9" and 9-12", cho- 
sen to have comparable area. The 1995 Keck data are 
of lower quality: the total integration was shorter, and at 
the time the instrument had much higher noise. Therefore, 
we used only the 3-9" annulus for this data. The model- 
subtracted galaxy image was multiplied by this mask, and 
the Fourier power spectrum of the central 256 x 256 pixel 
region determined. 

The power spectrum of the model-subtracted galaxy im- 
age has two components: (1) rising power at low wavenum- 
bers due to the fluctuations convolved by the point spread 
function (PSF), and (2) white noise due to the Poisson 
shot noise. The power spectrum P(k) can be represented 
as 



P(k)=P x£(k)+Pi, 



(1) 



0"5 FWHM seeing. Images on the galaxy were interlaced 
with blank sky fields and taken in an ABBA pattern. 

We obtained wider-field -ftf'-band (1.9-2.3 /zm) images 
on 28 June 1995 UT using the facility ncar-IR camera 
LIRC2 (Gilmore et al. 1995) on the Shane 3-m telescope 
at Lick Observatory. LIRC2 employs a Rockwell Interna- 
tional 256 x 256 HgCdTe NICMOS-3 array and has a plate 
scale of 0"38 pixel -1 . Data were obtained in a similar fash- 
ion as the Keck data, taking interlaced pairs of sky and 
galaxy images. Conditions were photometric with a typi- 
cal seeing FWHM of l'.'l. For all the runs we observed the 
faint IR standards of Casali & Hawarden (1992) as flux cal- 
ibrators. Hence, our resulting magnitudes are Vega-based. 

Details of the data reduction appear in Paper II. We sub- 
tracted an average bias from the images. Using twilight 
sky images, we constructed flat fields with an iterative- 
fitting algorithm which separated the flat field from the 
non-uniform thermal emission on the array. A prelimi- 
nary sky subtraction was performed to identify astronom- 
ical objects. Then for each image, we made a running sky 
frame from the prior and subsequent images of blank sky, 
excluding any astronomical objects. The sky brightness 
changes temporally so we scaled the running sky frames to 
the median counts in unsubtracted galaxy images, exclud- 
ing a circular region centered on the galaxy. We used the 
same masks for scaling the subtraction of the blank sky 
images. This ensured the blank sky images and galaxy 
images were reduced identically. NIRC images of bright 
sources suffer from "bleeding" : they have a positive hor- 
izontal trail which exponentially weakens along the read- 
out sequence. We used images of bright stars to model the 
bleeding and remove it. We also used a software mask to 
exclude a faint in-focus ghost of the exit pupil (the sec- 
ondary mirror) present in all the images. We used the 
galaxy to register the individual frames and averaged to 
assemble a final mosaic. 

Finally, we used the Lick images to determine the re- 
maining DC sky background in the Keck images. Being a 
cD-type galaxy, NGC 4874 has a rather shallow light pro- 
file; hence the common technique of assuming the galaxy 
follows a de Vaucouleurs r 1 / 4 profile is ill-suited. We ex- 
tracted azimuthally averaged profiles from the Lick and 
Keck images, and fitted for the DC level in the Keck im- 
ages. To estimate the offset between the slightly different 
Lick and Keck filters, we used spectra of solar-metallicity 
M0-M5 giants from Pickles (1998), which should be well- 
suited for this purpose (Frogel et al. 1978). Our synthe- 
sized (K' — K) color was —0.015 ± 0.003 mag, including 
accounting for the filters' fc-corrections at the redshift of 
Coma. 3 We checked our final calibrated images against 
the photometry of Persson et al. (1979) taken in a 14"9 
diameter aperture. The agreement was excellent, with our 
photometry being on average 0.014 ± 0.008 mag brighter. 

3. SBF MEASUREMENTS 

Our SBF measurement methods are similar to Tonry 
et al. (1990). A complete discussion appears in Paper II. 

The mean surface brightness of NGC 4874 was modeled 
by fitting for the harmonic content of isophotcs. Globular 

3 Our synthesized color is notably different than what would result from the transformation given in Wainscoat & Cowie (1992), which is 
K' — K = (0.18 ± 0.04)(_ff — K). However, as they point out, their transformation was derived from a sample of A stars and M dwarfs. Neither 
of these have the strong 2.3 (im CO absorption band seen in giant stars, which dominate the near-IR light of elliptical galaxies. The CO feature 
is within the standard X-band filter, but not in the _R"'-band filter; hence, the sign of our synthesized correction is as expected. 



where Pq is the total variance per pixel, Pi is white noise, 
and E(k) is the expectation power spectrum. P(k) is the 
convolution of the PSF power spectrum with the power 
spectrum of the software mask times the square- root of the 
galaxy model. It accounts for the PSF, the radial variation 
in the SBFs, and the effect of the software mask on P(k) 
(see Liu 2000). Using P(k), we fitted the two-dimensional 
power spectra to solve for P and Pi (Fig. 1). 

For the very lowest wavenumbers, flat-fielding and sky- 
subtraction errors produce extra power which contami- 
nates the SBF signal. We examined the power spectra 
of the blank sky fields, which were observed and reduced 
in an identical fashion to the galaxy images. We found 
that wavenumbers of k < 20 (0.08 pixel -1 ), corresponding 
to spatial scales of > 4x FWHM, had significant rising 
power. Therefore, we used only k — 20 to k = 128, the 
Nyquist frequency, when fitting the galaxy power spectra. 

We analyzed the power spectrum of each quadrant of 
the galaxy image independently We used the average of 
these four fits to obtain P and computed the standard 
error as the uncertainty in Pq. In Paper II, we ran Monte 
Carlo tests to verify that these quadrant-derived errors are 
accurate. Each Keck data set had only one suitable PSF 
star. Thus we adopted a PSF mismatch error of 0.08 mag, 
which we determined in Paper II from our Fornax cluster 
sample. 

Pq is the sum of the variance from the galaxy's stars, 
which is the desired signal, along with the variance from 
unresolved astronomical sources and instrumental effects. 
We quantified the effect of unresolved globular clusters 
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Fig. 1. — Tf-band fluctuation power spectra for NGC 4874. The galaxy power spectra were fitted by the sum (solid line) of a 
scaled version of the PSF (Po x E(k)) and a constant (Pi); the dashed lines show the contributions from these two components. 
We did the fit with the two-dimensional power spectra, using spatial frequencies of 0.08 pixel -1 to the Nyquist frequency. One- 
dimensional azimuthal averages are plotted to represent the results, with the measurement annuli and year of observation labeled. 



using an analytic representation for the GCLF, namely a 
Gaussian with three parameters (width, peak magnitude, 
and normalization) as measured by Kavelaars et al. (2000) 
using the HST U-band data. To convert to iGband, we 
adopted a color of V—K = 2.28, the average of the M31 and 
Milky Way globular clusters (Barmby et al. 2000). (The 
optical colors of the NGC 4874 clusters are similar to the 
Milky Way population, being relatively blue and having 
a moderately narrow dispersion [Harris et al. 2000].) We 
then computed the variance per pixel due to the globu- 
lar clusters below the detection limit of the HST optical 
images, (Jq C , using eqn. 5 of Blakeslee & Tonry (1995). 
Note that without the HST data, the variance from the 
undetected globular clusters would be comparable to that 
from the stellar SBFs, the signal we seek. 

Instrumental signatures from flat-fielding and sky sub- 
traction errors add variance to the lowest wavenumbers. 
By restricting our fits to k > 20, we reduced the system- 
atic errors. However, these effects might still contaminate 
the fitting, since they are not step functions in the power 
spectrum. This is particularly true for NIRC, which has 
spatially complex flat fields and bias frames. To quantify 
the amount of leakage into the wavenumbers used for fit- 
ting, wc used the same techniques on the blank sky fields 
and measured the blank sky variance per pixel o 2 k at 
numerous locations on the images. We used the scatter 
in the results to find the error in a 2 k , around 10-15%. 

While most of <J 2 sk arises from instrumental signatures, 
some comes from real astronomical sources, namely back- 
ground galaxies fainter than the detection limit. There- 
fore, unlike the usual SBF analysis, there was no need to 
compute the residual variance due to unresolved galaxies 
as this was included in (J 2 k . (Using the counts of Bershady 
et al. 1998, the variance from the galaxies is predicted to 
amount to «10% of the total variance.) 

Finally, the total amplitude of the power spectrum Pq 
was corrected by subtracting the contaminating variance 
P r . The quantity P r is the sum of Pqc an d P s ky, which are 
<7q C and a 2 k divided by the mean galaxy surface bright- 
ness per pixel in the measurement region, respectively. 
The exceptional depth of the HST optical data meant that 



residual variance from the unmasked globular clusters was 
negligible («0.002 mag). The blank sky variance was sub- 
stantial, about half of the total variance measured in the 
galaxy images. However, the measurements appear to be 
robust. For the 1998 Keck data, a 2 k measured in the 
two different annuli showed excellent agreement. In the 
1995 Keck data, a 2 k was substantially larger as expected 
by the higher noise levels then. After subtracting this 
variance, the resulting K-b&nd SBF apparent magnitude 
(WIk) agrees well with those from the 1998 data. 

The remaining quantity (Po — Pr) is the variance due 
solely to the stellar SBFs. To convert to mx , we used the 
measured photometric zero point and applied two correc- 
tions. We corrected for 0.003 mag of extinction (Schlcgel 
et al. 1998). Using the solar- metallicity models of Paper I, 
we also applied very small ^-corrections to tuk and V — I c 
to account for the redshifting of the galaxy light. The final 
errors on mx comprise the quadrature sum of the errors in 
the photometric calibration (0.02 mag), PSF uncertainty 
(0.08 mag), and measurement errors in the variances (oq C , 
a 2 k , and Po). Table 1 presents our results. 

4. RESULTS 

In order to compute the distance modulus, we use our 
ii-band SBF calibration from Paper II. This calibration 
uses 24 early-type galaxies in nearby clusters with high- 
quality SBF data; the zero point is based on /-band SBF 
distances to the bulges of six nearby luminous spiral galax- 
ies (Tonry et al. 2000) which have HST-derivcd Ccpheid 
distances from Ferrarese et al. (2000). (The effect of us- 
ing the new Cepheid distances by Freedman et al. 2001 is 
~0.05 mag — see Paper II.) The calibration is for the K$- 
band (2.0-2.3 /im; McLeod et al. 1995) filter, which has a 
slightly bluer bandpass than the K-bsmd filter used here. 
To estimate the effect of this difference, we use the theo- 
retical models of Paper I with ages of 3-12 Gyr and metal- 
licities of solar and slightly sub-solar (Z = 0.02 and 0.008); 
these span the SBF observations to date (Papers I and II) . 
The mean (Ks — K) offset is effectively zero (0.004 mag) 
with an rms of 0.02 mag. Thus, we assume no difference 
but include the rms as an additional source of error. The 
resulting calibration of the ii-band SBF absolute magni- 
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tude is: 

Mjf= (-5.84 ±0.04) + (3.6 ±0.8)[(V-J C ) - 1.15] . (2) 
This gives us the expected Mk for NGC 4874 given its (k- 
corrected) V — I c color. A weighted average of the distance 
moduli from the three data sets is 34.99 ± 0.15 mag. How- 
ever, the intrinsic (cosmic) scatter about the mean relation 
is non-zero; in Paper II, we conservatively estimated it to 
be 0.15 mag. We add this in quadrature to determine a 
final distance modulus of 34.99±0.21 mag (100±10 Mpc). 

Our errors do not include any systematic errors in the 
HST Cepheid distance scale. There are several poten- 
tial sources of error including the assumed distance to the 
Large Magellanic Cloud, the metallicity dependence of the 
Cepheid period-luminosity relation, and the photometric 
calibration of the HST WFPC2 instrument. Any changes 
in the HST Cepheid scale will affect our results by chang- 
ing the if-band SBF calibration. 

Our result is comparable to previous K-band SBF de- 
terminations of the distance to Coma. Thomsen et al. 
(1997) used 7-band observations from HST to derive a dis- 
tance of 102 ± 14 Mpc to the Coma elliptical NGC 4881. 
Jensen et al. (1999) measured a if-band SBF distance 
of 85 ± 10 Mpc to the other Coma supergiant galaxy 



NGC 4889. Both of these previous results were based on 
observations with much lower S/N, (P n — P r ) < 0.7Pi, at 
least a factor of 10 worse than our measurements here. 

To determine Hq, we adopt a mean velocity for the 
Coma cluster in the reference frame of the cosmic mi- 
crowave background of 7186 ± 428 km s _1 from Han & 
Mould (1992). At small redshifts, the H Q is related to 
the distance modulus by (m — M) = 25 + 51og(cz/Ho) + 
1.086(1 — qo)z. The last term accounts for cosmological 
curvature, and its effect on Hq is negligible: < 2% at the 
redshift of Coma for qg ranging from -1 to +1. Our re- 
sulting Hubble constant is 71 ± 8 km s _1 Mpc -1 . 
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Table 1 
NGC 4874 Ty-band SBF Measurements 



Date 


tint 

(s) 


Radius 

(") 


(mag/D") 


V-/ C a 
(mag) 


Pac/Po h 


Psky/Po 


(P -P r )/Pi c 


m K d 
(mag) 


18 Mar 1995 

19 Mar 1998 
19 Mar 1998 


1920 
3600 
3600 


3-9 
3-9 

9-12 


16.68 
16.68 
17.49 


1.219 ± 0.015 
1.219 ± 0.015 
1.198 ± 0.015 


0.00 ± 0.13 
0.00 ± 0.14 
0.00 ± 0.10 


0.69 ± 0.17 
0.48 ± 0.20 
0.58 ± 0.16 


6.7 ± 1.4 
9.3 ± 2.2 
7.2 ± 1.3 


29.52 ± 0.32 
29.47 ± 0.22 
29.18 ± 0.24 



a Mean K-b&nd surface brightness and V — I c color of the region used for SBF measurement. The V — I c data are from 
Harris et al. (2000) and have had fc-corrections applied. 

b The fractional contamination to the total SBF signal p). Pgc is the variance from globular clusters which are undetected 
in the HST optical imaging. P s k y is the variance measured from images of blank sky, which arises from instrumental noise 
and unresolved background galaxies. Both Pgc and P s k y are normalized to the mean surface brightness in the measurement 
region. The quoted errors include the error in Po, which dominates the error in Pgc / Po- 

c The signal-to-noise of the stellar SBF signal. P r is the total residual variance and is the sum of Pgc and P s ky The 
difference (Po — Pr) is the stellar SBF variance. Pi measures photon shot noise, the flat component of the power spectra plots. 

d The final SBF apparent magnitudes. The values have been corrected for extinction and have had fc-corrections applied. 
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